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Abstract. The Color Glass Condensate provides a systematic power counting of 
initial state effects in high energy QCD. We briefly discuss in this framework quark 
and gluon production in high energy collisions. 



1. Introduction 

The Color Glass Condensate (CGC) is an effective theory of the hadronic wave- 
function at high energies. For a review and references, see Ref. pQ. Key features, 
as summarized in the McLerran-Venugopalan model (MV- model) are as follows. 
The dynamical degrees of freedom are small x partons with occupation numbers of 
0(a~ 1 ) best described by the classical gauge fields A£. The large x partons are static 
light-cone sources of color charge represented by the two dimensional color charge 
density p a (x 1 _). The classical field of the nucleus obeys the Yang-Mills equations 
[D^, F^ u } a = 5 u+ 5(x~) p a (x±). The sources p a are described by the gauge invariant 
weight function Wa; [p a ] where xq separates "fields" and "sources". As xo decreases, 
fields are transformed into sources - the weight functional obeys a renormalization 
group equation (the JIMWLK equation) describing its evolution with x$. The evolution 
of the system is characterized by the saturation scale Q s (x) 3> Aq CD . One determines 
a posteriori that a s (Q s ) <C 1. Observables 0[A M ] are computed in the classical field 
for fixed p a and then averaged over W XQ [p a ] to obtain the gauge invariant expectation 
value: (O) = J[Dp a ]O[A^p a )]W X0 [p a ]. 

In the CGC framework, the ratio of the density p of hard parton color charges 
to the transverse momentum scale k± of interest (p/Ar[) enables a systematic power 
counting of initial and final state effects in the scattering of both dilute and dense 
partonic systems £Q. These will now be discussed. 



2. Gluon and Quark production in the dilute/pp regime (p p i/k\ , p p ilk\ <C 1) 

Inclusive gluon production in the CGC is computed by solving the Yang-Mills equations 
[D^, F^ u ] a = J%, where J v = 5(x~)5 u+ p p i + 5(x + )5 u ~ p P 2, with initial conditions given by 
the Yang-Mills fields of the two nuclei before the collision. To lowest order in p v \jk\ and 
Pp2/kj_, one can compute inclusive gluon production analytically. This was first done in 
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the A T = gauge [3] and subsequently in the Lorenz gauge d^A^ = 4j. The result 
at this order is fc^-factorizable into the product of the unintegrated gluon distributions 
in the two projectiles \. The comparison of this result to the collinear pQCD gg — > gg 
process and the fc^-factorized gg — > g was performed in ref. [Hj- This result for gluon 
production is substantially modified, as we shall discuss later, by high parton density 
effects in the projectiles. 

The result for inclusive quark pair production can be expressed in fc^-factorized 
form as |H| 



do~\ 1 f d 2 ki± d 2 k 



dy p dy q d 2 p ± d 2 q ± (2n)W A J (2tt) 2 (2tt 



2 | 



X if 1 (k 11 _)ip 2 (k 2 ±) ^-T2 , (!) 



where <z>i and (p 2 are the unintegrated gluon distributions in the projectile and target 
respectively (with the gluon distribution defined as xG(x,Q 2 ) = J d(k 2 J (p(x, k±)). 
The matrix element Tr k 2 ; q,p)\ 2 is identical to the result derived in the k±- 

factorization approach, which has been applied extensively to study heavy quark pro- 
duction at collider energies |7j|. In the limit ki± , k 2 ± — > 0, Tr \m~J~(ki, k 2 ; q,p)\ 2 /kf ± kl ± 
is well defined - after integration over the azimuthal angles in eq. d one obtains the 
usual matrix element \M.\ 2 gg ^ qq) recovering the lowest order pQCD collinear factorization 
result. 



3. Gluon and Quark production in the semi-dense/pA regime 

(p p /kl < i,p a A1~ i) 

Here one solves the Yang-Mills equations [_D M , F^ u ] = J u (with J v = 5 u+ b~(x~) p p (x±) + 
5 u ~5(x + ) p A (x±))) to determine the gauge field produced at lowest order in the proton 
source density and to all orders in the nuclear source density. The computations are 
performed in Lorenz/covariant gauge d^A^ = 0. Gluon production, in this framework, 
was first computed by Kovchegov and Mueller jS]. In ref. 0, the gluon field produced 
in pA collisions was computed explicitly. One obtains 

A"(q) = A^q) + J ^{Cfaku.) [U(k 2± ) - (2n) 2 5(k 2± )} 

+ C^q) [V(k 2± ) - (2n) 2 6(k 2± )]}^^ , (2) 

with k 2 = q — ki and U & V Wilson lines containing all orders in the nuclear source 
density p A . The coefficient functions and are simply related to the well known 
Lipatov effective vertex C M through the relation C M = C 1 + |C£. 

The path-ordered exponential U is a color matrix arising from the rotation of the 
color charge density of the proton source due to multiple scattering off the nucleus. The 

% This quantity is not the usual unintegrated distribution but is closely related. See Ref. JT] for a 
discussion. 
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path-ordered exponential V (differing from U by a factor 1/2 in the argument of the 
exponential) arises from the propagation of the produced gluon through the nucleus. 
Interestingly, the V^'s do not appear in the final result for gluon production. This is 
because for gluons produced on-shell one finds remarkably that C v ■ C v = C 2 = 
and = C 2 L = Akf^^/q^- Thus only bi-linears of the Wilson line U survive in 
the squared amplitude that gives the gluon production cross-section. The result is 
/c^-factorizable, except that now one replaces <£>2 with the unintegrated nuclear gluon 
distribution ip> A oc (WU). This distribution contains powers of the usual unintegrated 
gluon distribution to all orders - one recovers the usual unintegrated gluon distribution 
(if2 in eq. at large transverse momentum. 

Our result in the Lorenz gauge is exactly equivalent to that of Dumitru & McLerran 
in the A T = gauge ^Oj- The Cronin effect in proton-nucleus collisions has been studied 
by us in Ref. [Oj and several other authors previously and since and will not be discussed 
further here. 

Quark production can now be computed with the gauge field in eq. El [llj. The 
field is decomposed into the sum of 'regular' terms and 'singular' terms; the latter 
containing a factor S(x + ). The regular terms are the terms where a) a gluon from the 
proton interacts with the nucleus and produces a gg-pair after the collision, b) the gluon 
produces the pair which then scatters off the nucleus. Naively, these would appear to 
be the only possibilities in the high energy limit where the nucleus is strongly Lorentz 
contracted. However, in the Lorenz gauge, one has terms in the gauge field (the 'singular' 
terms, proportional to S(x + )) which correspond to the case where the quark pair is both 
produced and re-scatters inside the nucleus! Indeed, the contribution of this term to the 
amplitude cancels the contribution of the term proportional to the V's in the regular 
part. 

Our result for quark pair production, unlike gluon production, is not k±- 
factorizable. It can however still be written in fc^-factorized form as a product of 
the unintegrated gluon distribution in the proton times a sum of terms with three 
unintegrated distributions, <p 9 A ' 9 , <p 99,9 and ip 99,99 . These are respectively proportional 
to 2-point, 3-point and 4-point correlators of the Wilson lines we discussed previously. 
For instance, the distribution (p 99 ' 9 can be interpreted as the probability of having a 
qq pair in the amplitude and a gluon in the complex conjugate amplitude. For large 
transverse momenta or large mass pairs, the 3-point and 4-point distributions collapse 
to the unintegrated gluon distribution, and we recover the result for pair production 
(eq. HJ) in the dilute/pp limit. 

Single quark distributions are straightforwardly obtained. Here, the 4-point 
correlator in <p q ^' q9 collapses to a 2-point correlator ip 9 ^ 9 corresponding to a quark (or 
anti-quark) in the amplitude and complex conjugate amplitude. 

For Gaussian sources, as in the MV-model, these 2,3 and 4-point functions can be 
computed exactly as discussed in ref. [TTJ. Single quark distributions in the MV-model 
were recently computed by Tuchin [T2"] . 
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4. Gluon and Quark production in the dense/ A A regime 

This case is the relevant one for particle production in heavy ion collisions. It involves 
solving the Yang-Mills equations to all orders in the sources of both nuclei. This problem 
has not been solved analytically thus far - fc^-factorization breaks down completely here, 
even for gluon production. 

The problem has however been solved numerically . Non-perturbative formulae 
are derived relating (for collisions of identical nuclei) the saturation scale Q s in the 
nuclear wave-function to the energy and number distributions of gluons produced 
immediately after the collision (on a time scale ~ 1/Q S ). The gluon distribution is 
infra-red finite and is fit by a massive Bose-Einstein distribution with a "temperature" 
T ~ 0.47 Q s and m ~ 0.04 Q s for k± < 1.5 Q s . For k± > 1.5 Q s , it is fit by the tree level 
perturbative form pf ln(47rfcj_/ 'Q s ). 

The classical field description is valid as long as the occupation number / is greater 
than unity. As the system evolves, it becomes dilute and the classical description breaks 
down. Baier, Mueller, Schiff and Son [14J estimated this time to be 0(a s 3 ^ 2 Qj 1 ) 
In their "bottom-up" scenario, they suggest that inelastic 2^3 processes, though 
parametrically suppressed, are actually more efficient in driving the system from the 
classical stage towards thermalization. The analysis is valid for very small couplings 
and suggests that thermalization may take several fermis to achieve. In this light, the 
early thermalization required in some RHIC phenomenology appears puzzling. Arnold, 
Lenaghan and Moore have suggested that collective instabilities might drive the system 
faster towards equilibrium [To] . 

An interesting possibility is suggested by the following. As the classical field 
expands, one identifies a scale A(to) at an early time r which separates high momentum 
particles from low momentum classical fields. The high momentum particles scatter off 
the fields while the classical fields interact with each other and with the particles. With 
time, at an appropriate tl, one can define a new ("coarse graining") scale A(ti) at 
which one re-defines field and particle modes. We have developed an algorithm for 
a scalar field theory which implements this dynamical coarse graining while ensuring 
energy- momentum conservation ^B|- While promising, much work remains to extend 
this formalism to gauge theories. 
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